Producing cluster states in charge qubits and flux qubits 
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We propose a method to efficiently generate cluster states in charge qubits, both semiconducting 
and superconducting, as well as flux qubits. We show that highly-entangled cluster states can be 
realized by a 'one-touch' entanglement operation by tuning gate bias voltages for charge qubits. We 
also investigate the robustness of these cluster states for non-uniform qubits, which are unavoidable 
in solid-state systems. We flnd that quantum computation based on cluster states is a promising 
approach for solid-state qubits. 
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One-way quantum computing T], which is based on 
a series of one-qubit measurements starting from cluster 
states of a qubit array, is an intriguing alternative to the 
widely studied approach using unitary quantum gates. 
Here, the power of quantum mechanics, such as quan- 
tum parallelism and entanglement, is already stored in 
the initial cluster state. Cluster states are fixed highly- 
entangled states that involve all qubits and act as a uni- 
versal resource for quantum computing. 

Because of their unique importance, cluster states have 
been studied in a variety of physical systems. They have 
been extensively explored in optical quantum computers 
both theoretically 2] and experimentally 3]. By incor- 
porating cluster states, optical quantum computing can 
achieve substantially simpler operations 2] compared to 
the original linear optics quantum computing proposal 
0. Cluster states have also been studied in solid state 
qubits. In particular, processes of generating cluster 
states for single and encoded spin qubits have been pro- 
posed 0, m using the Heisenberg exchange interaction. 

The existing methods of generating cluster states all 
require multiple steps because of the types of interaction 
involved (in the case of photonic qubits, a large number 
of optical elements is also required). Here we describe 
theoretically an efficient method to create scalable clus- 
ter states in charge qubits 7, 8, 0,0,0, 0,0] and flux 
qubits using existing Ising-like interactions. 

Our key result is that cluster states in charge qubits can 
be created by applying a single gate bias pulse, right after 
preparing an initial product state |^o) = = 0)) — 
Hi|-h)i, where \±)i = (|0)i±|l)i)/V2. We also calculate 
the time-dependent fidelity of the cluster states in charge 
qubits using a quantum dot (QD) system with decoher- 
ence produced by the measurement back-action, and ex- 
plore the effects of non- uniformity among qubits, which 
is a realistic characteristics for all solid-state qubits. 

Cluster states in charge qubits. — The Hamiltonian for 
an array of charge qubits with nearest-neighbor interac- 



tions is described by 



Hcq — ^(f^j <^ix + i-l Oiz) + ^ Jlj CTi2 Oj 



(1) 
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with Pauli matrices and for the i-th qubit. fJ; 
is either the inter-QD tunnel coupling for coupled QD 
systems iSj iSi : or half the Josephson energy for super- 
conducting charge qubits [ll], [ij, |l3 , respectively. For 
either semiconducting or superconducting charge qubits, 
ti is the charging energy, and corresponds to the energy 
difference between |0) and |1) for each qubit. The cou- 
pling constants Jij are derived from the capacitance cou- 
plings. In one-way quantum computing 1], the Oix term 
needs to be switched off during the creation of the cluster 
state (Hcs), and then switched on when measurements 
are carried out. From this perspective, charge qubits 
with tunable fli 0| are desirable. However, tunability 
can produce decoherence and cross-talk between qubits 
themselves and between qubits and the environment. In 
addition, for most qubit systems, once the qubit array 
is made, fli and Jy are fixed, and only is controllable 
via the gate voltage bias (we call these "simple-design 
qubits"). Practically, such a simple design is preferable 
for solid-state qubits so as to simplify fabrication and 
enhance scalability. Thus, our goal here is to generate 
cluster states for simple-design qubits. Hereafter, with- 
out loss of generality, and for convenience, we focus on 
qubits of coupled QDs. |0) and |1) refer to the two single- 
qubit states in which the excess charge is localized in the 
upper and lower dot, respectively (see Ref. [qI). 

Typically, cluster states are generated by an Ising-like 
Hamiltonian Hcs = (ff/4) X]i<j(l^^iz)(l~o-jz)' where i 
and j are nearest-neighbor sites, starting from an ini- 
tial state \^o). Preparing a unitary evolution Ucs{t) = 
exp(itHcs) (we use h — 1) at gt — (2ni + 1)tt, where ni 
is an integer, is the first step for one-way quantum com- 
puting. Since simple-design charge qubits have Ising-like 
interactions, all we need to do to get Hcs out of Eq. 
is to turn off the effect of the a.^ terms in ^ . This can 
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be achieved in the high-bias regime where ^ fii > Jij , 
by applying a canonical transformation to Eq. [l^ : 

i i<j 

(2) 

where S = -« I]J^»/(2ei)]CTiy , with fl/{2ei) < 1. 
Also, Si = Ei + ^i/^i, neglecting terms of 
higher-order than (fli/ei)'^. H^w is an mi- 

wanted interaction term given by i?uw = 
~ YltLi Y.ia(.^i/'^i)Jhi+yd ^ix ^i+^d,z for a d-dimensional 
qubit array. Here 71 = ±1, 72 ~ {(±1, 0), (0, ±1)} and 
73 = {(±1,0,0), (0, ±1,0), (0,0, ±1)} are nearest neigh- 
bor indices for one-, two- and three-dimensional qubit 
arrays, respectively. As long as ifuw is sufficiently small 
and can be neglected, we can periodically generate 
cluster states in the tilted frame \^{t)) = e~'^|\E'(t)) 
after a time ics, if both Jt^s — tt/A + 2nj7r and 
Ei tcs = — (7r/4)ni -|- 271^7: are satisfied (n^ is the number 
of qubits connected to the i-th qubit; nj (> 0) and 
TT-E are arbitrary integers, and Jij should be uniform: 
Jij = J)- These two equalities lead to the relation 
J(8nE — ni)/{8nj + 1) = Ei = Ei^t^s- Thus, to generate 
a cluster state, gate bias voltage e,; for the i-th qubit 
needs to be set at 

e.^.f^%+^(%)'-nf m 

during tcs = 7r(8nj + 1)/(4J). To ensure a solution for 
e'i^ exists, we require tie — (2rii/J)nj > {2Qi/J + hi)/8. 
Hereafter we choose to treat the case of nj = 0, which 
corresponds to the shortest possible time to generate 
cluster states. To further ensure the high-bias regime, 
where 3> S^i, we require 8nE — ?ii ^ 2Qi/J. 

An initial product state |5'(0)) — |^'o) has to be pre- 
pared for qubits other than the input qubits P|. |^(0)) = 
H,(cos[(r!,/(2e,)+^/4] |0) +sin[(r!,/(2e,)+ ^4] |1)) is the 
corresponding state in the original {|0), |1)} basis, which 
can also be adjusted by the gate bias on each qubit. 

For example, to obtain the two-qubit cluster state 
I*)c3 = (|0)i|±)2±|l)i|-)2)/y2, where ni = fiz = 1, 
we can use Eq. (O and choose a gate bias (e^'' = e^**) 
from {6.37J, 14.7J, ...} {m > 1) for Q = 2J, and 
{13.8J, 22. 3J, ...} (nE > 2) for Q = 4J, etc. 

Our approach is valid as long as the unwanted 
term H^^ can be neglected. We can estimate a life- 
time, beyond which we lose the cluster state due to 
the presence of H^w, by calculating the fidelity de- 
fined by F{t) = (5'o|e^-f^-*e-'(-f^-+^™)*|*o) ~ 1 - 
zt(*o|ifuw|*o)+(l/2) (itf (*o|[ifuw,ifcs]+i?Ll*o)- For 
a d-dimensional A^-qubit array, 

F(t) « 1 - (^jJtj ^dN (4) 

Thus, the lifetime of the cluster state is limited by t < 
tuw = (2J(ri/e)\/d]v')^^. Furthermore, the constraint 



^cs < ^uw imposes a limit on the number of clustered 
qubits: A^^ax < {2e / {■k^))'^ / d. For example, consider a 
one-dimensional qubit-chain with 51^ = 4 J, and fii — 2. 
For riE = 2 [ef / J « 12.7), ^,„ax = 4. For ue = 6 
/ J ~ 45.6), A^max = 52. Indeed, we can choose an 
infinite number of bias conditions for each set of fixed 
and J. These are closely related to the possible number 
A^max of clustered qubits and to the scalability of the 
system. Various kinds of errors, as discussed in Ref. 0], 
should be taken into account for more detailed estimates. 

As shown in Ref. , Jij and Q.i are determined by the 
distances between QDs. Jy is directly determined by the 
capacitance network of QDs and basically controlled as 
a linear function of the distances between QDs. 17^ de- 
pends exponentially on the distances between two QDs 
in a qubit. With current technology it is quite difficult 
to fabricate an array of QDs with very uniform U,i. For 
superconducting charge qubits the situation is similar. 
However, notice that our approach does not require Sl^ 
to be extremely uniform since we can adjust according 
to Q.i in order to obtain an appropriate Ei in Eq. Q . In 
addition, as noted above, tie can be selected arbitrarily, 
which adds flexibility to our scheme. Thus, the one-touch 
cluster state generation method we discuss here should 
work with any charge qubit architecture. In short, al- 
though in general charge qubits have shorter decoher- 
ence times compared with spin qubits, our simpler and 
faster generation method could make them competitive 
with spin qubits in the context of cluster states, since 
several steps are required to generate cluster states for 
spin qubits 

Measurement scheme in charge qubits. — In one-way 
quantum computing, calculations are carried out by a 
series of local measurements in the Gx and Gy eigenbasis. 
For most charge qubits, however, the measurement is car- 
ried out in the cr^ eigenbasis {|0), |1)} by simply applying 
a large gate bias and using field-effect detectors such as 
quantum point contacts (QPCs) or single-electron tran- 
sistors. Thus, for charge qubits the ax and Uy basis mea- 
surements should be converted into cr^ measurements af- 
ter rotating the frame via ±(7r/2)j^ and ±(7r/2)a; pulses. 
These pulses can be generated by applying AC gate bi- 
ases such as ei{t) = eoiCos{LUict + (pi). In a rotating co- 
ordinate ?7rw(0 = exp(— i Wic^cTi^:), the wave function 
is given by |5'(t)) ~ U^^Y^{t)), and the Hamiltonian on 
resonance {ujic — ^i) in given by: 

i?rw ~ ^ ^('^'^ coS(^i - aiy sin(j)i) + Hy^ , (5) 

i 

where Hy^^J^iji^v /'2)i'^iy '^jy+'^tz cTjz) is an unwanted 
term here. In order to realize |^'(tm)) — e^'-^^i"'^ l*(0)) 
at the measurement time t = tm, we should have 
gi//rwt„, ^ e*S'"'*»'^"eT^Ti:>'^'!/. Thus, we obtain <„! = 
{TT/il)li, the voltage amplitude e^o — (7r/2t„i)(l ± 4/2), 
and the phase 4>i = 7r/2 -I- l^n {li, I2 and ^3 are arbitrary 
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integers {h ^ 0)) for the Ux measurement. If we take 
<m = (7r/(4f)))(l + 4/1) and Cio = 27r?2/im (^2 + 0), we 
realize the Oy measurement. By applying a gate bias af- 
ter the time tm, and Oy measurements are processed. 
This scheme works well for the region where ti^j J 1. 

For the experiment in Ref. ^^Ii ^ = 20.5 /icV and J = 
95 /ieV, and thus ^cs — ^/(4^) = 34.2 psec, ~ 664 /icV 
from Eq. © at tie = 1 and t^-^ = (tt/ J)(ecs/[2V2ri]) = 
1.57 nscc for a dephasing time of the order of T2 ~ 5 
nsec. Thus, our approach can be applied, if smaller J 
(< f2) is prepared, for instance, by increasing the distance 
between qubits. 

Cluster states in flux qubits. — For flux qubits, fl > e 
in Eq. (Q) , so that here we cannot directly use the above- 
mentioned one-touch approach for charge qubits. Here 
we show a method to generate cluster states for flux 
qubits by applying an oscillating magnetic field. Con- 
sider two inductance-coupled flux qubits working at the 
optimal bias |l6l| with Hamiltonian: 



Hfq = ei (Tl^ + £2 0-22 + f^f COs(Wi t + (jjl) (Tla 
+ COs(w^'i + 02) Cr2x + J'xx CTlx <^2x , 



(6) 



where ilf^ and Lof are the half amplitude and the fre- 
quency respectively of the applied classical field. At 
the optimal bias point, the system is immune, up to 
first order, to variations on the control parameters and 
is thus robust against decohcrence. This Hamiltonian 
is a good starting point for generating cluster states. 
In the rotating wave approximation for two identical 
qubits (r^i = ^-i), we have iJfq ~ Hq + H^y, with 
^0 = Z]Li(^f cos0i -f fj^j^ sin0j) and H^y = 
Jxx{pixO'2x + o'iyO'2y)- The Operator to generate cluster 
states, Uca, is produced by switching on and off the reso- 
nant field of fif- and controlling the phase (f>i similarly to 
the conventional conditional phase gate operation. For 
example, if we define Ria{9) = exp(z0CTiQ,), (a = x,y), 
and Uxy{0) = exp(z6'((Tia;tT2a -I- uiyU2y))^ we have: 

Ucs = Rlx{()l)R2x{Q2)Uxy{0z)Rlx{04^Uxy(0^Rlx{0Q) 

= exp(-i7r((Ti^+0-22; -0-ij;0-2a;)/4) (7) 

with Oi = 62 ^ -7r/4, 6*3 = 6I5 = 7r/8, 6*4 = 7r/2 and 
= — 7r/2. After rotating Ucs around the y-axis, we re- 
cover the original cluster state generator Ucs- In the case 
of many qubits, cluster states for the entire system are 
generated by applying Eq. Q to all the neighboring qubit 
pairs. Note that Eq. ((HJ also describes the rotating wave 
approximation to the charge qubits in Eq. Thus, this 
method of generating cluster states is also applicable to 
charge qubits, and not only to flux qubits. 

The time required for the creation of a cluster state in 
flux qubits is Tflu^ = 57r/(2r2f ) + 7r/(4f]^) + tt/{Uxx)- 
Taking fjf - 17^ ~ J^, - 0.5 GHz, we obtain Tflu^ ^ 18 
nsec (T2 ~ 200 nsec [iSj). The effect of imperfect pulses 
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FIG. 1: Time-dependent fidelities F(t) of cluster states 
|*>C2, |*>C3, |*)c4 , a product state |1010) and four-qubit 
decoherence-free (DF) states |*)df4 = (|1100) - |1001) - 
jOllO} + |0011»/2, for To = J and AF = 0.6J. f7 = 4J 
and riE = 2 in Eq. @, thus £2 = £3 ~ 12. 7J, ei = £4 ~ 13. 8J. 
(a) Comparison of the cluster states, DF state and |1010). (b) 
The case when non-uniformity in the qubit parameters is in- 
troduced as VLi = 4J(l-r;i), u = e+rjiJ and F,'=^' = (l-77i)r<*' , 
with i indicating the i-th qubit. Here 77^ = for all qubits 
besides rjs = 0.1 (i) and 774 = 0.1 (ii). The fidelities of |^')c4 
for (i) and (ii) mostly overlap, (inset) Four qubits that use 
double dot charged states are capacitively coupled to a QPC 
detector. We consider a similar detection setup for two and 
three qubits. These calculations are carried out by the Hcq, 
that is, they include Huw and higher order terms. 



If we take the deviation from a perfect pulse as Si — S2, 
63 — S^, 64 = 5e, for the initial state |^'o), the fidelity is 
given by 

F{t) w l-5f/2-2Sl + sm^[n/8] {Sl/2) + iSi6ie'''/^ (8) 

Thus, the fidehty F{t) remains close to one, up to second 
order in the pulse shape error, even when the pulse shape 
has defects. 

Effect of non-uniformity in cluster states. — Cluster 
states are highly entangled states involving all qubits, 
and can be decohered by various kinds of local fluctu- 
ations. Here we investigate the effect of non- uniform 
qubit parameters on cluster states in semiconductor QDs, 
from Eq. using a measurement setup, which pro- 
duces decoherence in the double dot qubits through back- 
action. We analyze a capacitively coupled detector (such 
as a QPC), whose shot-noise constitutes a random charge 
fluctuation on the qubits. 

We use a density matrix (DM) to describe up to four 
qubits (inset of Fig^l The DM equations for the 

qubits and the QPC detector are derived in Ref. [l9| : 



(9) 



pl/2_pl/2 

Zl Z2 



can be estimated by substituting 6j 



+ Sj in Eq. |7J). 



where zi,Z2 = (1111), (1110), (0000) for four qubits 
(256 equations) and zi,Z2 = (11), (10), (01), (00) for two 
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qubits (16 equations). K^^ is the energy of the Zi state 
and depends on and Jy in Eq. For example, for 
two qubits, if(ii) =ei+e2 + >^i2 while i^(io) =£1^ e2^>^i2- 
The gj{zi)s are introduced for notational convenience 
and are determined by the relative positions between 
qubit states. We assume that the tunneling rate F of 
the QPC detector in the presence of N qubits satisfies 
F"-'^ = , where the tunneling rate F^ is deter- 

mined by the state = ± 1 of the i-th qubit. The 
strength of the measurement can be parametrized by 
ATi as f|^^ = F,o ± AF„ where F io is the tunneling 
rate of the QPC in the absence of the qubits. The 
time-dependent fidelity F{t) = Tr[/i(0)/5(t)] can be cal- 
culated by tracing over the elements of the reduced DM 
obtained from Eq. I^. F{t) can be expanded in time 
as F{t) = l-E,rf(l/"-0(*/'^*"'')"> where the lifetime is 
^ {_Tr[p(0)d"p(0)/cir]}i/". 
From Eq. 0, we obtain the first-order lifetime 
for two-, three- and four-qubit cluster states |^')c2 7 
I*)C3 = a+)i|0)2|+)3 + |-)i|l)2|-)3)/v^, and = 

(|+)l|0)2| + )3|0)4 + | + )l|0)2|-)3|l)4 + |-)l|l)2|-)3|0)4 + 

|— )i|l)2|+)3|l)4)/2, respectively, as follows: 

l/r« ^ ^d(zi,.2)/8 (10) 

Zl,Z2 = (ll),..,(00) 

^l^cl = E A(^i,z2)/32, (11) 

zi,Z2 = (lll),..,(000) 

1/^4 = E rd(zi,Z2)/128, (12) 

zi,Z2 = (llll),..,(0000) 

1/2 

where the dephasing rate is defined as / d(zi, 22) = [Fzj — 
Fz2^]^. Note that the lifetime of the cluster states is an 
average over all the dephasing rates between different 
product states. This is in contrast with other entan- 
gled states. For example, the lifetime of two-qubit Bell 
states |c) = (|10) + |01))/\/2 and \d) = (|10) - |01))/\/2 
takes the form 1/ri^^ = l/r^^ = (1/2)^(10,01). It is 
well known that the singlet state jd) is the most robust 
two-qubit state |2nl | when there is a symmetry between 
qubits. However, solid-state qubits generally decohere 
due to various kinds of local causes, which often break 
the symmetry of the qubit state. Our results in Eq. (10- 
12) indicate that cluster states might be robust against 
non-uniformity or local defects because of the averages. 

In Fig. n] we compare the fidelity of cluster states with 
a product state and a four-qubit decoherence-free (DF) 
state 20]. Figure n^a) shows the time-dependent fideli- 
ties of two-qubit and four-qubit cluster states and the 
product state |1010), when Fq = J and AF = 0.6 J. Our 
results show that the strongly entangled cluster states are 
more fragile than a product state such as |1010). We can 
also see that the robustness of the cluster state depends 
on the number of qubits in the cluster state. Figure^b) 
shows the time-dependent fidelities of both cluster states 



and DF states, for non-uniform qubits. Here the param- 
eters ti and F^ for the third or fourth qubit deviate 
from those of other qubits by 10%. Note that the fideli- 
ties of |^')c4 show almost the same behavior irrespective 
of the distribution of the non-uniformity. Furthermore, 
a comparison between Fig.^a) and Qlb) shows that the 
non-uniformity has almost no effect on the fidelity of the 
|'I')c4 cluster state. These results vividly illustrate our 
analytical analysis of the lifetime (that it is an average 
over all the product states). In contrast, in Ref. we 
showed that the robustness of the DF states strongly de- 
pends on the non- uniformity. Thus, even though cluster 
states are generally more fragile than DF states, they are 
more robust against non-uniformities among qubits than 
DF states. 

In conclusion, we describe how to efRciently gener- 
ate cluster states in solid-state qubits. By manipulating 
the gate bias voltage, we explicitly show how to gener- 
ate 'one-touch' entanglement via cluster states in charge 
qubits. We also investigate the robustness of cluster 
states, and find that one-way quantum computing could 
be viable for solid-state qubits. 

FN and XH are supported in part by NSA, LPS, ARO, 
and NSF. 
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